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We study analytically the time evolution in decaying 
chaotic systems and discuss in detail the hierarchy of charac- 
teristic time scales that appeared in the quasiclassical region. 
There exist two quantum time scales: the Heisenberg time tn 
and the time tq = tn / x/kT (with k ^ 1 and T being the de- 
gree of resonance overlapping and the transmission coefficient 
respectively) associated with the decay. If tq < tn the quan- 
tum deviation from the classical decay law starts at the time 
tq and are due to the openness of the system. Under the oppo- 
site condition quantum effects in intrinsic evolution begin to 
influence the decay at the time tH- In this case we establish 
the connection between quantities which describe the time 
evolution in an open system and their closed counterparts. 
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In a recent paper [Q G. Casati et at, using numerical 
simulations in the kicked rotor model with relaxation, 
have demonstrated that a new time scale exists for a de- 
caying quantum system in the deep quasiclassical region. 
After this time, which is much less than the Heisenberg 
time tn = 'inp (with p being the mean level density and 
h = 1), the decay law begins to deviate from the classical 
one. 

The aim of the present paper is to show that such a 
time scale is, in fact, a general feature of open quantum 
chaotic systems and is related to peculiarities in fluctua- 
tions of the resonance widths. We describe these fluctu- 
ations in the framework of the random matrix approach 
and employ the formalism of the effective non-Hermitian 
Hamiltonian which is commonly used in the theory of 
resonance scattering [^,^. 

Generally, decay properties of open quantum systems 
are related to fluctuations in complex eigenvalues (reso- 
nance energies) = En — (i/2)r„ of the effective Hamil- 
tonian Ti. via the two-point correlator of the Green's op- 
erator g{£) = {£ — 'H)~^ . As typical examples one can 
mention the ^-matrix m § or time delay §-0,|ll cor- 
relation functions. The simplest quantity of such a kind 
is the leakage of the norm inside an open system 



P{t) = (1/iV) (Tr{exp(iHtt)exp(-i7^t)}) 



(1) 



which is a somewhat simplified version of the decay func- 
tions considered in Ref. Q. Here, the angle brackets 
stand for the random matrix ensemble average, and the 
equivalence of spectral and ensemble averages is implied 
| p2[ . It is easy to see that similar quantity (without the 
ensemble averaging) has been numerically evaluated in 



§. In the case of a closed system P{t) = 1 at any time. 
The time dependence in eq. (^ appears due to antiher- 
mitian part of the operator Ti. The well-known relation 
between the time evolution operator exp(— iTii) and the 
Green's function Q{E) enables one to represent P{t) in 
the form of the Fourier integral 



oo oo 



-oo — oo 



In the eigenbasis of the effective Hamiltonian the trace 
in Eq. (^ can be represented in the form 



= ^\E^«'" exp{-z(£„-£:,)0) (3) 
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where _R(e,r) denotes 



= ^ (E ^n'J [e - {En ~ )\5[T- 



and Un'n = (V'n'IV'n) is the Bell-Steinberger nonorthogo- 
nality matrix of the eigenvectors oi Ti. This 
matrix differs from the unity only if resonances over- 
lap. It is worth noting that, contrary to the e depen- 
dence of the function i?(e, F) determined by the level 
spacings along the real energy axis, its F dependence 
is governed by widths themselves. Therefore, the de- 
cay law can not be directly related to the distances 

{En — En'Y + (F„ — F„')2/4 between resonance levels 
in the complex energy plane. This is contrary to what 
was conjectured in Ref. 

Prior to exact calculating Pit) we would like first to 
perform qualitative analysis. As it will be justified below 
by the exact calculation, the main features of P{t) can 
be understood already from calculation of the diagonal 
part 



Pd{t) 



OC 



Tt 
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Here we put approximately = 1, neglecting its 

smooth dependence on the index n. The function ViT) 
is the distribution of the resonance widths. This hmc- 
tion is explicitly known |^,^ for the case of the unitary 
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ensemble which corresponds to the systems with the bro- 
ken time-reversal symmetry. Therefore, we use this en- 
semble to demonstrate our general statements. In Ref. 
I psf detailed analysis of the width distribution has been 
performed, in particular, convenient for our purposes in- 
tegral representation 
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Krf {M-iy. 
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(6) 



),k(1-T)/T 



is given, with rj = T/Tw being the decay width measured 
in the units of the Weisskopf width [ p^JT?) 



w 



MT/2TTP , 



(7) 



where T is the transmission coefficient and the dimcn- 
sionless parameter n = 2iTpTw = MT characterizes the 
degree of resonance overlapping. The rate of small widths 
diminishes rapidly when the number M of (statistically 
equivalent) open decay channels grows. For small over- 
lapping, K ^ 1, the density V{r]) simplifies to the well- 
known Xai distribution. However, quasiclassics corre- 
sponds to M ^ 1 and strong overlapping k ^ 1 . In 
this case the width distribution decreases exponentially 
at r < Tw and follows the power law ~ (Tw /r)^ within 
some domain above Tw ■ In the classical limit M, p —^ oo 
but Tw is kept fixed and identified with the classical es- 
cape rate 0, an empty strip appears below the value 
Tw Hi- 

Substituting Eq. (|) in Eq. (|) and changing the order 
of integration, we come to the expression 



T/(l-T) 
i f 

T J (iHY 



exp 



-Mln[l 
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Twt 



(8) 



which is still exact in M and T. In the classical limit 
defined above the first term in the 1/M expansion of 
the logarithm in Eq. (||) gives Pd{t) — Pci(t)p{t), where 
Pci{t) ~ exp{~Twt) is the classical decay probability 
which follows from the semiclassical periodic orbit theory 
| p9[ and p(t) is a slowly varying factor, the proper calcu- 
lation of which lies beyond the diagonal approximation. 
Further terms of the 1 /M expansion can be neglected for 
the times appreciably less than the characteristic time 



'Mtw — \J n/Ttw = tH/vKT 



(9) 



where tw = ^/Tw is the characteristic lifetime of the 
system. The quantum time scale tq is similar to that 
found in ref. j^] (see also |^). In the mesoscopic sys- 
tems the typical life time is given by the Thouless time 
pl[ . Therefore, the connection tH = Ktw shows that 
our overlapping parameter k plays the role analogous to 
the dimensionless conductance in the mesoscopic physics. 
We note in this respect that the ratio tw /tq differs from 



that conjectured in Ref. |^ by an additional factor \/T, 
which depends on the strength of coupling to channels. 
It is also worth noting that the time tq appears in the 
relaxation phenomena in disordered conductors as well 
|2|,||. 

The next-to-leading term of the expansion being pos- 
itive, quantum corrections slow down the decay law at 

t > tq 

totic power law 



oq. After this time crossover occurs to the asymp- 



p(-)(t) = «-i (Twt/M) 



-M 



(10) 



which is characteristic for open quantum systems pJ^ . As 
it will be demonstrated below, this expression correctly 
matches the exact result. The fact that the diagonal ap- 
proximation properly reproduces the asymptotic behav- 
ior was first noted in |Q . One can easily see from Eq. ^ 
that such a power behavior comes from the infiuence of 
the widths which are smaller than Tw- Their rate differs 
from zero as long as the parameter 1/M remains finite. 

Qualitative arguments presented above can be put on 
a rigorous ground. Powerful supersymmetry technique 
p5|,p| enables us to perform exact calculations in Eq. (j|). 
Skipping the details of quite a standard calculation, we 
concentrate on the analysis of the result. For the case of 
unitary symmetry it reads 



Pit) = dXo dXi p{K) /(AO 5{t/tH - (Ai-Ao)/2) 



l + r(Ao-l)/2 



l+T(Ai-l)/2 



M 



(11) 



where p,{\i) — (Ai — Ao)~^ is the measure of integration 
pSf and f{Xi) = (A^ — Ao)/2. The openness of the sys- 
tem is contained in the last "channel" factor in Eq. (|Tl]). 
Actually, the structure of the expression (|l^) is of univer- 
sal nature for different quantities that describe the time 
evolution of a chaotic quantum system. It consists of 
the integration with the measure which is specific for the 
chosen ensemble, the channel factor, and the preexponent 
f{Xi). The latter is the only factor that depends on the 
concrete quantity considered. Since the time dependence 
related to decay properties comes just from the channel 
factor, our analysis is of quite a general meaning. 

At t < tff the decay probability (11) can be repre- 
sented in the form 



1 

P{t) = jdv[l + {l- 2u)t/tH\ 
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whereas dX t > t^ it looks like 



(12) 



2 



P{t) 



Mlnil+r^t/M) /di. [1 + (1 - 2u)tH/t] 
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iyT/il + Twt/M)\ ■ ^^^^ 

In the classical limit (see above) simple calculation leads 
to the classical decay law Pci{t). 

The peculiarities of quantum deviation from the clas- 
sical time evolution depend significantly on the ratio 



'kT. If kT > 1, so that < < tn, 



the analysis of Eq. dlj) goes along the similar way as de- 
scribed below Eq. (|P^and leads to the same conclusions. 
In the opposite case kT <C 1 (but still k ^ 1, which im- 
plies small values of the transmission coefficient T) one 
arrives, by inspecting the integral factor in Eq. (p^, at 
the general relation 



Popenit) = {l+Twt/M)" 



PrI 



lit). 



(14) 



which remains valid till the time tf — tw/T — tq/y/ nT — 
tn/nT. Here Pciosed{t) is the Fourier transform of the 
spectral correlation function which describes the time 
evolution in the corresponding closed system. Analogous 
calculation of the function P{t) for the case of the or- 
thogonal ensemble (which corresponds to chaotic systems 
with time-reversal symmetry) gives instead of Eq. ( [Tl| ) an 
expression of similar structure but with obvious changes 
which are characteristic to the symmetry class consid- 
ered. We only mention that the channel factor contains 
the power M/2 rather than M . Therefore, the same re- 
lation ([l^) with M substituted by M/2 is valid also for 
time-reversal invariant systems. 

The decay factor (1 -t- Vwt/M)~^'^ is equivalent to the 
classical exponent till the time tq. However, in the taken 
case the Heisenberg time tu is smaller than tq, which 
yields the influence of quantum effects on the time evo- 
lution via the function Pdosedit)- At last, after the 
time Ktf — tn/T the asymptotic regime P^°'^\t) ~ 
{Twt/M)~^''^ appears with a proportionality coefficient 
depending on the quantity considered. In the case of the 
function (^ its asymptotics coincides with that given by 
Eq. (0). 

In conclusion, even in the quasiclassical domain there 
exists a finite probability of the widths less than the clas- 
sical escape rate Tw- This leads to the appearance of the 

Ih/V^ asso- 



new quantum time scale tq = n/Ttw 
ciated with the decay. The parameter of resonance over- 
lapping K,':$> 1 plays the role analogous to the dimension- 
less conductance in condensed matter physics. The quan- 
tum effects begin to influence the time evolution starting 
from the time tq if tq/t^ = I/a/kT ^ 1 and from the 
Heisenberg time tu under the opposite condition. In the 
latter case the relation (^^ holds connecting the time 
evolution in an open system with its closed counterpart. 
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